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ABSTRACT 

We study the signifieance of major-merger-driven star formation in the early Universe, 
by quantifying the contribution of this process to the total star formation budget in 80 
massive (M* > 10^° M©) galaxies at z ~ 2. Employing visually-classified morphologies 
from rest-frame F-band HST imaging, we find that 55^^''^% of the star formation 
budget is hosted by non-interacting late-types, with 27^®% in major mergers and 
18^^% in spheroids. Given that a system undergoing a major merger continues to 
experience star formation driven by other processes at this epoch (e.g. cold accretion, 
minor mergers), ^^27% is an upper limit to the major-merger contribution to star 
formation activity at this epoch. The ratio of the average specific star formation rate 
in major mergers to that in the non-interacting late-types is ~2.2:1, suggesting that 
the enhancement of star formation due to major merging is typically modest, and that 
just under half the star formation in systems experiencing major mergers is unrelated 
to the merger itself. Taking this into account, we estimate that the actual major- 
merger contribution to the star formation budget may be as low as ^^15%. While our 
study does not preclude a major-merger-dominated era in the very early Universe, if 
the major-merger contribution to star formation does not evolve strongly into larger 
look-back times, then this process has a relatively insignificant role in driving stellar 
mass assembly over cosmic time. 

Key words: galaxies: formation - galaxies: evolution - galaxies: high-redshift - 
galaxies: star formation - galaxies: interactions - galaxies: bulges 



1 INTRODUCTION 

Understanding the processes that build massive galaxies is 
a central topic in observational cosmology. The observed 
peak in the cosmic star formation rate (SFR ) at z ~ 2 (e.g. 
iMadau et allll998l : [Hopkins fc Beacomll2006l ') indicates that 
a significant fraction of the stellar mass in today's massive 
galaxies is likely to have formed around this epoch. However, 
the principal mechanisms that created this stellar mass re- 
main unclear. Was the star formation driven by vigorous, 
major-merger (mass ratios > 1:3) induced starbursts? Or 
were processes other than major mergers - e.g. cold accre- 
tion, minor mergers, etc. - responsible for creating the bulk 
of the stars in today's mas sive galaxies, as suggested by 
recent theoretical work (e.g. iKeres et aflboogl : lOekel et al.l 
l2009l )? 

Modern surveys that access large UV/optically-selected 
samples of galaxies at 2: > 1.5 have facilitated the empiri- 



cal study of star formation around z ~ 2 (e.g. [Paddi et al.l 
2004 lErb et al.ll2006l:[Reddv et al.ll2005l: lOaddi et al ll2007l : 



Santini et all l2009l : iHathi et all I2OI0I : IWuyts et all 12011 
Star-forming galaxie s at this epoch lie on a star-formatio n 
'main sequence' (e.g. lOaddi et al.ll2007l : iReddv et al.ll2012l '). 

where galaxy SFRs are proportional to their stellar masses 
with a slope of unity (relatively passive galaxies lie be- 
low this sequence). The growing body of observational 
work on these galaxies increasingly suggests that much 
of the cosmic star formation at this epoch may be un- 
related to the major-merger process. Integral-field spec- 
troscopy of star-forming galaxies around 2 ~ 2 has re- 
vealed a high fraction of systems with properties indica- 
tive of turbulent disks and only a modest inc i dence of ma- 
or m e rgers (e.g. Forster Schreiber et al. 20061: Genzel et al] 



20081: IShapiro et 8,1.1 l2008l: iForster Schreiber et al. 



Cresci et al.ll2009l : ICenzel et al.ll201ll : iMancini et all 



200S; 



20111 . 
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Figure 1. Example J+H composite images for the morphological 
classes into which we split our galaxy sample. Galaxies are classi- 
fied into spheroids (top row), major mergers (disturbed systems 
with multiple nuclei and clear, extended tidal features, middle 
row) and non-interacting late- types (bottom row). An angular 
scale which spans half the image is shown. The images are ~45 
kpc on a side. Note that the images usually appear better on 
screen than in print. 



see also Law et al. 2009, van Dokkum et al. 2011). Imaging 
studies, that have explored the rest-frame UV and optical 
morpholog i es of s tar-forming galaxies at t hese e pochs (e.g 



es 01 s tar-tormmg galaxies at t hese e pocns (e.g. 
I2OO6I : iForster Schreiber et all I2OIII : iLaw et al.l 
^K^..7^o indicated a preponderanct^ non-mergmg 



Lotz et al. 



systems amongst high-redshift star formers, suggesting that 
the role of major mergers may indeed be subordinate to that 
of other processes in d riving star format i on at this epoch. 

In a recent study, iRodighiero et all (|201ll ) have shown 
that 'starbursts' - systems that show enhanced star forma- 
tion and lie off the main sequence of normal star-forming 
galaxies - have a relatively minor role at this epoch, ac- 
counting for around 10% of the cosmic star formation ac- 
tivity. However, starbursts can be driven either via ma- 
jor mergers or by dense nuclear star-forming regions (e.j 
iDi Matteo et al.H2007l : iDekel et al.]|2009l : iDaddi et al1l201C ' 
More importantly, many major mergers share the same star- 
formation chara cteristics as norm al star-forming galaxies at 
this epoch (e.g. iLaw et al]|2012al . see also Di Matteo et al. 
2007, Kaviraj et al. 2012) and thus lie on the main star- 
formation sequence itself As a result, a unique one-to-one 
mapping is unlikely to exist between major mergers and star- 
bursts. To probe the relative significance of major-merger- 
driven star formation at z ~ 2, it is desirable to quantify 
the proportion of the total star formation budget that is 
attributable to systems that are morphologically selected as 



major mergers at this epoch. This has not been directly ad- 
dressed by previous work and represents both a quantitative 
empirical result and a useful constraint on theoretical mod- 
els at high redshift. 

Deep near-infrared imaging from current WFC3 sur- 
veys - which trace rest-frame optical wavelengths at 2 ~ 2 
- enables us to morphologically classify massive galaxies at 
this epoch and study how star formation activity is appor- 
tioned in terms of galaxy morphology (e.g. major mergers, 
non- interacting late- types, etc). It is worth noting, however, 
that a system undergoing a major merger at 2 ~ 2 contin- 
ues to experience star formation driven by gas inflow via 
other processes such as cold accretion and minor mergers 
(major mergers can be thought of as simply the 'clumpiest' 
part of the mate rial flowing in alon g the cosmic web). Both 
theoretical (e . g. iDekel et al.l [2OO9I) an d observational (e.g. 
iKavirai etlo] I2OI2I : iLaw et all l2012al ) work indicates that 
star formation due to these other processes is significant at 
this epoch and possibly comparable to the major-merger- 
driven activity. Hence, in addition to splitting the star for- 
mation budget by morphology, it is necessary to consider 
the proportion of star formation in major-merging systems 
that is unrelated to the major-merger process, and subtract 
this from the fraction of the star formation budget hosted 
by systems with major-merger morphologjQ. 

Here, we probe these questions using a complete, rest- 
frame optically-selected sample of massive (M. > 10'° M0) 
galaxies at 2: ~ 2, drawn from the WFC3 Early Release 
Science (ERS) programme, which provides unprecedentedly 
deep near-infrared HST imaging and ten-filter photometry 
in the GOODS-South field. Section 2 describes the galaxy 
sample that underpins this study. In Section 3, we describe 
the derivation of galaxy properties e.g. SFRs, stellar masses 
and internal extinctions. We study the proportional contri- 
bution of major mergers to the total star formation bud- 
get in Section 4 and summarise our findings in Section 5. 
Throughout, we use t he WMAP7 cosmological parameters 
iKomatsu et al.ll201ll) and pres e nt ph otometry in the AB 
magnitude svstem loke fc GunnI (|l983l ). 

2 GALAXY SAMPLE AND 

MORPHOLOGICAL CLASSIFICATIONS 

The WFC3 ERS programme has imaged ~45 arcmin^ of the 
GOODS-South field in the WFC3 UVIS (F225W, F275W, 
F336W) and IR (F098M [Y], F125W [J], F160W [H]) chan- 
nels, with exposure times of 1-2 orbits per filter. The ob- 
servations, data red uction, and instrument p erformance are 
described in detail in IWindhorst et all (l201l[). Togeth er with 
the existing ACS BViz imaging ( Giavalisco et al.] 20041 ). the 
data provide 10-band panchromatic coverage over 0.2 - 1.7 
/^m, with 5(T point source depths of AB <^ 26.1 — 26.4 mag 
in the UV and AB ^ 27.2 - 27.5 mag in the IR. 

Here, we focus on an _ff-band selected sample of 80 ERS 
galaxies, that have stellar masses M, > 10'° M0 and red- 
shifts in the range 1.9 < z < 2.1. For 12% if our sample 
(10 galaxies) we use published spect r oscop ic redshifts from 
ISantini et all (|2009l ): IPopesso et aD (|2009h : iFerreras etaP 



^ Note that the situation is significantly different at low redshift, 
where gas-rich major mergers can enhance star formation by or- 
ders of magnitude fe.g. lMihos fc Hernauistll996l ). because secular 
processes drive star formation weakly. Almost all the star forma- 
tion in low-redshift major mergers is, therefore, attributable to 
the merger itself. 
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l|2009l ). For the remaining objects we use photometric red- 
shifts, calculated using the EAZY code (Brammer et al. 
2008). The peak of the redshift probability distribution from 
EAZY is used as the best estimate of the redshift - the ac- 
curacy of EAZY redshifts at this epoch is Az ~ 0.1 and 
the nominal time interval our defined by redshift range is 
^0.3 Gyr. The H-hand traces rest-frame V at 2: ~ 2 and the 
galaxy sampl e is complete within the se stellar mass and red- 
shift ranges (|Windhorst et al]|201ll ). Studying the massive 
end of the galaxy population restricts us to systems that are 
both bright {H{AB) < 24.2 mag) and extended, which facil- 
itates reliable morphological classification. The narrow red- 
shift interval minimises both morphological K-corrections 
and overlap between the spheroid and major-merger mor- 
phological classes (as we discuss in Section 4). 

Here we classify galaxies via visual inspection of their 
composite J-KH images. Since the J and H filters corre- 
spond to the rest-frame optical wavelengths at 2; ~ 2, 
these images trace the underlying stellar populations in 
each galaxy and not just the UV-emitting star-forming 
regions. Visual classification of morphologies in the high- 
redshift Universe has been commonly employed in the lit- 
erature, using rest-frame optical HST images that have 



the ERS images used here (e.g. 


Windhorst et al. 20021: 


Cassata et al.l 2010l; Kavirai et al. 




2011: Cameron et al.l 


201ll;lKocevski et al.ll2012l: iLaw et al 


l2012al). Visual classi- 



fication offers better precision and consistency than morpho- 
logi cal parameters (such as CAS, M20, Gini coefficient, see 



e.g. lAbraham et al.ill996l: 
l2004r iTavlor-Mager et al.l 



Conselice et al.ll2003l: iLotz et all 



20071 ) . which can be mo re sen- 



sitive to image resolution a nd signal-to-noise (e.g. iLisked 
l2008l : iKartaltepe et al] I2OI0I ) but are valuable for classify- 
ing large datasets, where visual classification is prohibitively 
time-consuming. 

Galaxies are classified into the following three broad 
morphological classes: [1] spheroids [2] non-interacting late- 
types and [3] major mergers, which are disturbed systems 
that exhibit multiple nuclei and clear, extended tidal fea- 
tures. The number fractions in classes [1], [2] and [3] are 
19%, 43% and 38% respectively. Figure 1 presents examples 
of objects drawn from each morphological class. 

3 PARAMETER ESTIMATION: STELLAR 
MASSES AND INTRINSIC STAR 
FORMATION RATES 

A variety of methods have been employed in the literature 
to derive galaxy SFRs. Calibrations (often based on sam- 
ples of local galaxies) can be used to convert X-ray, UV, 
infr ared or radio luminosities into estimates of SFR (see 
e.g.lReddv fc Steidel200i : lDaddi et al.l2007l : |Pannella et a .1 



I2OO9I : Elbaz et al.l 201 ll ). Alternatively, galaxy spectral 
energy distributions (SEDs) can be fitted to theoreti- 
cal star formation histories (SFHs) t o derive SFRs (e.g. 
IShaplev et"al] l2005l : iLaw et al.l l2012al 'l. Since dust is in- 
cluded as a free parameter in these SFHs, intrinsic (i.e. 
dust-corrected) SFRs can be derived self-consistently using 
this method. The derivation of reliable SFRs ideally requires 
rest-frame UV photometry (as is the case here), since the 
leverage in the SFR and the dust extinction comes largely 
from these wavelengths. 

In this paper, we calculate galaxy SFRs via SED fitting. 
The WFC3/ACS photometry of each individual galaxy is 
compared to a large library of synthetic photometry, con- 
structed using constant SFHs, each described by a stellar 
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Figure 2. TOP: SFR as a function of galaxy stellar mass in 
our galaxy sample. The main sequence of star-forming galax- 
ies at 2 ~ 2 is shown using the solid line (the observed scat- 
ter is indicated using the dotted lines). Galaxy morphologies are 
shown colour-coded (black = spheroids, blue = non-interacting 
late- types, red = major mergers). BOTTOM: The distribution of 
derived extinction values for our galaxies. The median value is 
shown by the vertical red line and the median from the recent lit- 
erature (see Tresse et al. 2007 and Cucciati et al. 2012) is shown 
by the blue dotted line. 



mass (M), age (T), metallicity (Z) and internal extinction 
(Eb-v)- We vary T between 0.05 Gyrs and the look-back 
time to jz = 20 in the rest-frame of the galaxy, Z between 
0.1 Zq and 2.5 Z© and Eb-v between and 1 mag. Syn- 
thetic magnitudes are gene rated by folding the model SFHs 
with the stellar models of iBruzual fc CharlotI (120031 ) with 
dust attenuation applied following ICalzetti et al.l (|2000^ . 
The likelihood of each model, exp(— x^/2), is calculated us- 
ing the value of x^, computed in the standard way. Estimates 
for parameters such as stellar mass, internal extinction and 
SFR are derived by marginalising each parameter from the 
joint probability distribution, to extract its one-dimensional 
probability density function (PDF). We use the median of 
this PDF as the best estimate of the parameter in question, 
with the 25 and 75 percentile values (which enclose 50% 
of the probability) yielding an associated uncertainty. Since 
dust is explicitly taken into account in this process, we de- 
rive intrinsic SFRs, free of internal reddening, directly from 
the SED fitting process. The derived internal extinctions, 
SFRs and stellar masses are uncertain by ^0.1 mag, ^0.1 
dex and ~0.2 dex respectively. 

In the top panel of Figure (2] we plot our derived SFRs 
vs. galaxy stellar mass. The SFR values for our star-forming 
galaxies are consistent with the star formation main se- 
quence at these epochs defined by the recent literature. 
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Galaxies that lie below this sequence are typically spheroids, 
which are relatively passively-evolving systems. Recall that, 
unlike studies that specifically target star-forming systems, 
the mass-complete sample employed here is not biased 
against galaxies with low star formation rates. In the bot- 
tom panel of Figure [21 we present the distribution of derived 
internal Eb-v values for our galaxies. The spread in our val- 
ues (0 < Eb-v < 0.5 mag) agrees w ell with that found by 
other studies (e.g. iLaw et a]ll2012a ) and the median of our 
distribution {Eb-v ~ 0.25 mag i.e. Apuv ~ 2.2 assuming 
Calzetti et al. 20 00) is in good agreement with the litera ture 
a.tz~2 (see e.g.lTresse et al.l (|2CI07l ), [Law et all (|2012bl ) and 
ICucciati et"aLl (|2012l . see their Figure 4)). 



4 THE MAJOR-MERGER CONTRIBUTION 
TO THE STAR FORMATION BUDGET 

We begin by exploring how star formation activity is ap- 
portioned in terms of galaxy morphology, by summing the 
derived SFRs of galaxies in each morphological class and 
considering the fractional contribution of these classes to 
the total star formation budget (Figure 3). We find that 
55^^*% of the star formation activity takes place in non- 
interacting late- types, with 27**% in major mergers and the 
rest (18*^%) in systems that have spheroidal morphology. It 
is worth noting that the proportion of star formation driven 
by morphologically selected major mergers calculated here is 
hig her than the corr e spond ing value derived for starbursts 
by iRodighiero et al" As we noted in the introduc- 

tion, this is due to the fact that many major mergers exhibit 
similar or only modestly-enhanced SFRs compared to nor- 
mal star-forming galaxies, lie on or close to the star-forming 
main sequence (see Figure [2] above) and are, therefore, not 
part of the more extreme starburst population. 

The predominance of non-interacting late-types in the 
total star formation budget indicates that major mergers 
are not the dominant mechanism driving star formation in 
massive galaxies at 2 ~ 2. Furthermore, as we noted in the 
introduction, systems undergoing major mergers continue to 
experience star formation via other processes (e.g. cold flows 
and minor mergers). Hence, 27% represents an upper limit 
to the major-merger contribution to the star formation bud- 
get. To improve our estimate, we consider the enhancement 
of star formation due to major merging, since this better 
represents the portion of the star formation activity that is 
directly attributable to this process. While measuring this 
enhancement is not possible in individual major mergers, we 
can estimate a typical value for the population as whole by 
comparing the mean specific SFR in the major mergers to 
that in the non-interacting late-types. 

The ratio of the mean specific SFRs in these two mor- 
phological classes is ~2.2:1 (major mergers : non-interacting 
late- types), implying that, on average, around half the star 
formation in major-mergers is likely driven by other pro- 
cesses. This relatively modest enhancement in star formation 
activity due to major merging is consistent with the findings 
of recent theoretical work (e.g. ICenll20l"ll ) and also empiri- 
cal studies that do not find significant differences between 
the SFRs of galaxies that are morpholo gically disturbed and 
those that are n ot at this epoch (e.g. Kavirai et al.l [20121 : 
iLaw et al]|2012al ). Thus, if around half the star formation 
in major mergers is unrelated to the merger itself, then the 
major-merger contribution to the total star formation bud- 
get is likely to be as low as ~15% (i.e. 27% x 1.2/2.2). 

Before we conclude this section, we briefiy discuss the 
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Figure 3. TOP: The fraction of the total star formation budget 
in massive (M» > lO^'' Mq) systems at z ~ 2 that is hosted by 
various morphological types. Major mergers account for less than a 
third (27%) of the total star formation budget, while non- interacting 
late-type galaxies host more than half of the star formation activity. 
Given that systems undergoing major mergers continue to experience 
star formation driven by other processes, 27% is an upper limit to 
the major merger contribution to the star formation budget. The 
actual contribution is likely to be ~15% of the total star formation 
budget (see text in Section 4 for details). BOTTOM: A pie-chart 
visuaUsation of the star formation budget apportioned in terms of 
galaxy morphology (Sph = Spheroids, LTG = Non-interacting late- 
type galaxies, Maj. Mgr = Major mergers). 

spheroid population in the context of the major mergers. 
We note first that the time interval spanned by our study 
(~0.3 Gyr) is shorter than the effective times cales (0.5-2 
Gyr) over which major me rgers coalesce (see e.g. iLotz et al] 
|2oo4 In ewman et al]l2012l ). so that the morphological classes 
do not overlap with each other. More importantly, however, 
iKavirai et al.l (|2012l ) have used high-resolution cosmologi- 
cal simulations to demonstrate that spheroids at 1 < 2 < 3 
that are remnants of recent major mergers (i.e. ones that co- 
alesced within the last ~0.5 Gyr) will exhibit clear tidal fea- 
tures at the depth of the ERS images. This study has further 
demonstrated that many newborn spheroids in this redshift 
range do not carry such morphological disturbances, indicat- 
ing that a significant fraction of these systems are not built 
via major mergers (in agreement with the results of recent 
theoretical work, e.g. Dekel et al. 2009). Around 15% of the 
spheroids in our sample show morphological disturbances 
and these galaxies account for ~3% of the total star forma- 
tion budget. While the spheroid and major merger classes 
do not overlap (as discussed above), it is clear that, even if 
we added the disturbed spheroids to the major merger por- 
tion of the star formation budget, our conclusions would re- 
main unchanged. Our analysis therefore indicates that major 
mergers contribute a relatively insignificant fraction (~15%) 
of the total star formation budget in massive galaxies at 
z ~ 2 and are not the principal driver of cosmic star forma- 
tion at this epoch. 

5 SUMMARY 

We have explored the significance of major mergers in driv- 
ing star formation at high redshift, by quantifying the con- 
tribution of this process to the total star formation budget 
in a sample of 80 massive (M» > 10^" Mq) galaxies at z ~ 2. 
We have found that 55*^*% of the total star formation activ- 
ity in massive galaxies at this epoch is hosted by late-type 
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galaxies that are not interacting with other systems, with 
27^*% in major mergers, and the rest in spheroids. 

Since systems undergoing major mergers continue to ex- 
perience star formation driven by other processes (e.g. cold 
flows and minor mergers), ~27% is an upper limit to the 
contribution of major mergers to the total star formation 
budget. To improve our estimate, we have considered the 
typical enhancement of star formation induced by a major 
merger, since this is the portion of star formation activity 
which is directly attributable to this process. We have esti- 
mated this enhancement using the ratio of the mean specific 
star formation rate in the major merger population to that 
in the non-interacting late- type galaxies. In agreement with 
recent observational and theoretical work, we have found a 
relatively modest enhancement (x2.2), which implies that, 
on average, just under half the star formation activity in ma- 
jor mergers is unrelated to the merger itself. This reduces 
the contribution of major mergers to the total star forma- 
tion budget in massive galaxies at this epoch to ~15%. Our 
analysis therefore indicates that the contribution of major 
mergers to the total star formation budget in massive galax- 
ies at 2 ~ 2 is relatively insignificant and this process is not 
the principal driver of cosmic star formation at this epoch. 

Since our study is based on instantaneous star forma- 
tion rates, it provides only a snapshot of the star forma- 
tion budget in massive galaxies at z ~ 2. Thus, while our 
data cannot rule out a merger-dominated era in the very 
early Universe, if the major-merger contribution to stellar 
mass assembly does not evolve strongly into earlier look- 
back times {z > 2), then major mergers are unlikely to be 
significant contributors to the overall buildup of stellar mass 
in the Universe. In forthcoming papers we will use morpho- 
logical analyses of large datasets such as CANDEL S - e.g. 
via projects such as Galaxy Zoo (|Lintott et al.|[2008l ). which 
uses 450,000-1- members of the public to visually classify 
large survey datasets - to comprehensively study the con- 
tribution of major mergers to the star formation budget as 
a function of stellar mass, environment and redshift. 
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